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Abstract 

•4— > . 

A possible explanation of the apparent disrepency between the theo- 
retical prediction and experimental measurement of the ratio of lifetimes 
T(Ab)/T(Bd) is that "spectator effects", which appear at 0(1 /ml) in the 
r^ • heavy quark expansion, contribute significantly. We investigate this possi- 

£> . bility by computing the corresponding operator matrix elements in a lattice 

simulation. We find that spectator effects are indeed significant, but do 
5_i ■ not appear to be sufficiently large to account for the full discrepency. We 



stress, however, that this is an exploratory study, and it is important to 
check our conclusions on a larger lattice and using a larger sample of gluon 
configurations. 

PACS: 12.15.-y, 12.38.Gc, 12.39.Hg 



1 Introduction 

At leading order in the heavy-quark expansion the decay rate of the heavy quark 
is independent of its parent hadron. In this letter we present the results of an 
exploratory study, in which we attempt to gain some understanding of the striking 
discrepancy between the experimental result for the ratio of lifetimes |l | 

T(Afc) -0.78 ±0.04 (1) 



T 



(B d ) 



and the theoretical prediction [fj], based on the heavy-quark expansion @-[0] 

0.98 + 0(1 /ml) . (2) 



r ( A ^ _nnv nn /,,,m 



In particular, we compute the contributions to the 0(1/ ml) term on the right- 
hand side of eq. (0) which come from "spectator effects", i.e. from decays in 
which two quark or antiquark constituents of the beauty hadron participate in 
the weak decay. These effects may be larger than estimates based purely on 
power counting would indicate as a result of the enhancement of the phase space 
for 2 — > 2 body reactions relative to 1 — > 3 body decays [FJ. 

Denoting the 0(1 /ml) contribution from spectator effects by A and following 
the analysis of ref. we find that [] 



A = -0.1735! + 0.195e 2 + 0.030L! - 0.252L 2 (5) 



where 



£ iK) = T23 K~ = -0.01 ±0.03, (6) 

j B m B <tm B 

, , 8 (B g \bLt a qqRt a b\B cl ) 

have been computed, using lattice simulations, in ref. ||. The two new variables 
introduced here are 

8 (A b \h»LqqrLb\A b ) 

L ^ m b) = TI3 K~ (8) 

J B m B ^rn Ab 

8 (A b \h^Lt a qq^Lt a b\A b ) 

LaK) = ISHTs 2^ • (9) 

Heavy-quark symmetry implies that there are only two matrix elements which 
need to be considered for the A&, in contrast to the four for S-mesons 0. 

We stress that this is an exploratory study. It is the first calculation of the 
matrix elements L 12 and provides a preliminary indication of whether spectator 
effects can reconcile the experimental measurements and theoretical predictions 

1 Equation (||) can be derived from eq. (39) of ref. B by the substitution 

B = -6Li (3) 

2 The coefficients of two of the matrix elements for B-mesons are so small that we don't 
include them in eq. (|g). 



for the ratios of lifetimes in eqs. (|1|) and (0). We perform the calculations with 
a static 6-quark and two (rather large) values of the mass of the light-quark 
and do not attempt to extrapolate the results to the chiral limit. Our results 
indicate that spectator effects are not negligible, although they do not appear to 
be sufficiently large to account fully for the discrepency. Specifically we find: 

/ -0.31(3) for am n = 0.74(4) 
Ll[mb) ~ { -0.22(4) for am^ = 0.52(3) , (iUj 

/ 0.23(2) for am n = 0.74(4) 
L 2 {m b ) - | Q17(2) for am ^ = 052(3) 5 I 11 ) 

with a -1 ~ 1.1 GeV. The corresponding results for the ratio of lifetimes are 
r(A 6 ) _/ 0.91(1) for am, = 0.74(4) 



r(B d ) ' \ 0.93(1) for am n = 0.52(3) . 

2 Perturbative matching 

In this section we briefly discuss the matching factors which are required to obtain 
the matrix elements of the continuum four-quark operators renormalised at a scale 
\i from those of the bare lattice operators computed in lattice simulations at a 
cut-off a -1 . The details of the calculation are presented in ref. ||. Here we simply 
summarise the main points required for the evaluation of the matrix elements L\ 
and L 2 in the MS scheme. 

We start by using the renormalisation group to relate the matrix elements Li 
and L 2 , defined in the MS scheme at two different renormalisation scales, \x = m b 
and fi = a^ 1 . Since the Wilson coefficient functions in the OPE expansion (|5|) 
have been evaluated at tree level only ||, we keep just the leading logarithms in 
the evolution equations so that P, 110] 



Lx(m b )\ [ 1 + 2^* 



L2(m b ] 



c F s 




(13) 



where 




9/2/3 

- 1 = 0.40 ± 0.04 . (14) 



In estimating 5 we have used A QCD = 250 MeV, a -1 = 1.10 GeV, m b = 4.5 GeV 
and (3q — 9. The error in 5 is evaluated includes a 20% uncertainty for Aqcd- 

In the second step of the matching we relate the matrix elements renormalised 
in the continuum to those regularized on lattice, both at the same scale, a -1 . 
Although this involves corrections of 0(a s ), which are, in principle, beyond the 



precision which we require, we nevertheless include them because the perturbative 
coefficients in lattice perturbation theory are generally large. The computation 
for the most general four quark operator involving one heavy quark appears in 
the appendix of ref. ||. For L\ and L 2 the relevant relations are: 

L ifi~ l ) = 72 [h u M 1 + h 12 M 2 + h 13 M 3 + h^M^} (15) 



P B m B 
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L 2 {a~ 1 ) = — 2 [/121M1 + h 22 M 2 + h 23 M 3 + h 2A M A ] , (16) 



where 



Mi = (A 6 1 (brLq) (jfLb) \A b ) 



2m k 



M2 = (Atl (bYLfq) fil" Lb) |A„) | (A 6 | (bj^Lq) (g 7 VL6) \A b ) 
2m A6 2m Ab 

Mg = (A 6 | gr^g) (qYLt a b) \A b ) 

2m Ab 
M ___ (A b \{h' 1 Li rq){qj> i Lt a b)\A b ) | (A 6 | (&7^Lt°g)(g7VLt a 6) |A 6 ) 

2m Ab 2m Ab 

are the matrix elements of the bare lattice operators regularised at a -1 . The 
coefficients hy are listed in table ([I]) where for the lattice coupling constant we 
have used a boosted coupling equal to 

a^ia' 1 ) 6(8/t crit ) 4 , , 

-^-T ~ = 71 To j - 0-01216 . 21 

47T (47r) 2 /? 

Readers who prefer to use other choices of the lattice coupling constant can 
combine the coefficients in table p] with their choice of coupling. 

A consequence of the Heavy Quark Effective Theory, and the fact that the 
light quarks in the A b are in a spin zero combination, is that the number of 
lattice operators whose matrix elements have to be evaluated is four rather than 
8 (which is the case for heavy mesons ||). 

In order to obtain the factor f B m B it is also necessary to determine the 
normalization of the axial current, Z A , 

f B m B = (0| Aoia- 1 ) \B) = Z A (0| ^(a" 1 ) \B) = JmT b Z a Z l , (22) 

where A l Q lt is the time component of the axial current defined on the lattice and 
Zl is obtained from the matrix element determined in numerical simulations At 
our value of the lattice spacing 

Z A = 1- 20.0^ '- ~ 0.75 . (23) 

In the following section we combine the results for the matrix elements computed 
on the lattice (M» and Zl) with the perturbative coefficients presented in this 
section (hij and Z A ), to obtain the values of L\(m b ) and L 2 (m b ). 



coeff. 


expression 


value 


/in 


1 + t If ~ 5*1 - !^J * 1 - 21-65^ 


0.737 


hu 


t [-1*3] * 9.19f 


0.112 


hi3 


^[-|-x 4 -2x 5 -x 7 ] ^9.29^ 


0.113 


hu 


£ [-x 6 ] =< 6.89f 


0.084 


tl21 


^|[-|-x 4 -2a; 5 -a: 7 ]^2.06f 


0.025 


h 2 2 


# [-1*6] * l-53f 


0.019 


h 2 3 


1 + 5 [A - l*i + 1*2 - Ix 4 + |x 5 - Jx 7 ] ^ 1 - 10.82£ 


0.869 


h 2 4 


£ [W ~ -3.45^ 


-0.042 



Table 1: Matching coefficients for the matrix elements M\ and M 2 , renormalised 
on lattice at an energy scale a -1 = 1.10 GeV. The values of the integrals Xi are 
reported in the Appendices of ref . M . 



3 Lattice computation and results 

The non-perturbative strong interaction effects in spectator contributions to in- 



clusive decays are contained in the matrix elements M{ in eqs. (17)— (2C). They 
have been evaluated in a quenched simulation on a 12 3 x 24 lattice at /3 — 5.7 
using the SW improved action |T| , 



S 



sw 



S 



gauge 



+ s 



Wilson 



icsw 



Y^ q(x)F flu (x)q(x) 



(24) 



X,fl,U 



where S sange and S Wllson are the Wilson gauge and quark actions, respectively. We 
use 20 gauge-field configurations and the light quark propagators are computed 
using a stochastic inversion based on the exact relation 



I* 



".1 



— j[d(f>]{Aj h (j) k )*(f)i exp (-$(A^A) lm <f) m ) 



(25) 



where the are auxiliary bosonic fields, introduced in order to perform the in- 
version of the matrix A, which in our case is the fermionic matrix. To reduce the 



statistical noise the technique of maximal variance reduction has been used [|12|] . 
The use of this stochastic inversion technique makes it possible to compute a 
light-quark propagator from each point in the half of the lattice with < t < 12 
(we call this region box I) to each point in the other half where 12 < t < 24 (box 
II). This increases considerably the effective statistics in the computation of the 
matrix elements. 

We have performed the calculations at csw — 1-57, which is the numerical 
value of 1/wjj, with w being the the average value of a link variable as defined from 
the trace of the plaquette. The calculation is therefore "tadpole-improved" and 
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Figure 1: Plot of Qf after the substraction for the contribution of excited states. 



hence the perturbative coefficients in table (|1|) are the same as those which would 
be obtained with a tree-level improved action (csw — !)• We have evaluated 
the matrix elements with two values of the light-quark mass, corresponding to 
Ki = 0.13847 (for which m^a ~ 0.74(4)) and n 2 = 0.14077 (for which m^a ~ 
0.52(3)) [0. For this value of csw K crit — 0.14351 |T^|. The same lattice has 
been used to compute, with satisfactory results, the wave function of a I?-meson 
and the effective coupling constant for the decay B* — ► B + it (in the Heavy 
Meson Chiral Lagrangian). This computation is reported in ref. |14| and the 
values which we use for ]\vhb are extracted from this paper. 

The evaluation of the matrix elements requires the computation of two- and 
three-point correlation functions of the form, 



C 2 (t x ) = J2(0\J(x)J\0)\0) 



(26) 



where we have assumed L. > 0, and 



c 3 (o,t x ,t y ) = J2(u\Ay)o(o)f(x)\o) 



j-,y 



(27) 



where t y > > t x . In eqs. (|26| ) and (p7[) J and J^ are interpolating operators 
which can destroy or create the A& baryon, for which we take 



j\=e abc (u a a {l*C) a /,) I 
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(28) 



where u, d, b are the quark fields. In eg. ([28]) a, b, c are colour labels, a,/3, 7 are 
spinor labels and a sum over repeated indices is implied. We define Z\ by 



Z A u { °\0) 



(A b ,s|Jt(0) 7 |0) 
\/2m A 



(29) 



where s represents the spin state (up or down) of the baryon. 

In order to enhance the contribution of the ground state to the correlation 
functions, it is useful to "smear" the interpolating operators J and J^ . In this 
paper we will follow ref. |12] and adopt the type of smearing known as "fuzzing" . 
This technique consists in replacing light quark field q(x), by a "fuzzed" field 



q F (x)= J2 U l F (x)q(x + i) + U 



Ii{x)q{x 



(30) 



i=l,2,3 



where U^x) are defined by the recursive relations 

U[(X) = P SU (2) 



U F Ax) 



V, 



SU(2) 



(U t F (x) + Y,Uf(x)U t F (x + ])U F 3 (x + i + j) + 

u^{x)uf(x-j)u[{x + i-: 

C U F t (x) + J^UfixW^x+JW^ - i + j) + 



U F 3 {x)U F t {x-])U F {x-\ 



(31) 



(32) 



starting with initial values U F {x) = U{(x) and U^x) = Uj(x — i). Vsu{2) is a 
projector on SU(2), implemented as in the Cabibbo-Marinari cooling algorithm, 
and ( = 2.5 is a constant value. The recursive procedure for U F t (x) has been 
applied twice. 

We introduce two superscripts on each correlation function, each of which can 
be either U F" or "L" , which indicate whether the interpolating operators J and 
J^ are fuzzed or local. 

The standard technique to extract hadronic matrix elements of the type 
(A&| O I Aft) is to look for plateaus in the ratios 



R(0,ti,h 



72 C( F (Q,t u t 2 ) 
' K Ck F {t x )C% F {t 2 ) 



(33) 



In our analysis, however, even with the use of fuzzed interpolating operators J 
and J\ we cannot eliminate the effects of excited states from the three-point 
correlation functions C[ F (0,ti,t 2 ) satisfactorily. On the other hand we do find 
that the ground state dominates the two-point correlation function C FF (t) for 
t > 3, and the masses we obtain in this way agree, within errors, with those found 
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Figure 2: Plot of the matrix elements Li(o x ) and L 2 (a x ) computed after the 
substraction of excited states. 



previously on the same lattice [12[ using a 3-mass correlated fit for a number of 
smeared correlators. In order to obtain the matrix elements M, of eqs. (|I"7D-([20|) 
we therefore need to subtract the effects of the excited states. 

We have followed the following procedure to extract the matrix elements 
(A b \0\A b ): 

• For each value of the light-quark mass we start by fitting the two-point 
correlation function for t > 3 with a single exponential 

C FF (t) = (Z F i ) 2 exp(-m A t), (34) 

thus obtaining the mass of the ground state, m A . Within errors, the 
masses of the ground state which we obtain are in agreement with those 
obtained from the same lattice using a more sophisticated fitting procedure 



in ref. 12 . 



• We model the contribution of the excited states by a second exponential 
and now fit C FF (t) for t > 1 by 

Ci F (t) = {Zife-^ + (Z£) 2 e— V , (35) 

keeping m A and Z F fixed at the values obtained from the single exponential 
fit above. We find that C FF is well represented by the two exponentials. 

• For each operator O we then fit the three-point correlation function 
C[ F (0,t u t 2 ) to 

C FF (0,t 1: t 2 ) = iM^M(^)V mA (l tl '+' 2 ) (36) 

2mA 

I Q I" — m^\tx\-rn\ x t2 i -m\\ti\— mA 1 *2] (?>7\ 

obtaining values for the two unknown parameters (A&| O |Aj,) and the con- 
stant C, which encodes the contribution from excited states. 

This procedure has been repeated for each of the 4 relevant operators, and for 
the linear combinations corresponding to L\ and L 2 on 40 jackknife samples to 
extract the statistical errors. 

In order to control the contributions from the excited states more effectively it 
will be necessary to carry out a simulation with considerably improved statistics. 
It is, however, possible to check the consistency of our approach a posteriori. We 
subtract the contributions from the excited states obtained above from the two- 
and three-point correlation functions, and look for plateaus in the ratios: 



expression 



Kl 



«2 



Z|Z| = /im B /2 


0.33(1; 


0.33(1) 


Mi 


-0.026(3; 


-0.019(3) 


M 2 


0.045(4; 


0.039(5) 


M 3 


0.018(2; 


0.013(2) 


M 4 


-0.040(4; 


-0.031(4) 


L^a- 1 ) 


-0.18(2; 


-0.13(3) 


L 2 (a- 1 ) 


0.21(2; 


0.16(2) 


Li{m b ) 


-0.31(3; 


-0.22(4) 


L 2 {m h ) 


0.23(2; 


0.17(2) 



Table 2: Lattice results for the matrix elements computed on lattice, Mi, the 



combined matrix elements at two different scales, 
lifetimes. 



Li, and the physical ratio of 



where the tilde indicates that contribution from excited states has been sub- 
tracted from the correlation function. The subtracted two-point correlation func- 
tion is reported in fig. [IJ Fig|J shows the plateaus for the ratios R corresponding 
to the operators in L\ and L 2 (with the appropriate normalization factor fi 8 m ). 
The plateaus in figfj give us confidence in our treatment of the subtraction of the 
contribution of the excited states. The results for the matrix elements obtained 
from eqs. ([J7]) and ( j38|) agree to within 1%. 

Our results for the matrix elements at each of the two values of k are reported 
in table 0. Combining them with eq. (0) we obtain 



r(A fe ) 
t(B°) 



0.91(1) 
0.93(1) 



for arrin 
for am n 



0.74(4) 
0.52(3) 



(39) 



^From eq. (|3"9"| ) we see that, although they are of 0(l/m\), spectator effects 
are indeed significant (compare eqs. (|39|) and (fj)). Estimates of the parameter r 
defined in eq. (HI), using the non-relativistic quark model or the bag model |IBL IT? 



or QCD Sum Rules |18|] are typically in the range 0.1-0.5. On the other hand, 
Rosner has estimated r from the spin splitting between Sq and Eq baryons 
(Q=c,b) and finds r ~ 1 (2) from charmed (beauty) baryons |l!|. A recent 
reanalysis of this problem using QCD Sum Rules in which more condensates are 
introduced, finds a range of possible values for the ratio of lifetimes (including 
ones close to the experimental value in eq.( |])), depending on the (unknown) 
values of various condensates f20(| . At the two values of the light quark mass 
at which we do our computations we find r ~ 1.2 ± 0.2, which is at the high 
end of expectations. The large values which we find for r and consequently the 
significant effect on the prediction for the lifetime ratios, make it important to 
improve the precision of the lattice simulations. 



10 



In quark models and related pictures, the parameter B = 1. In our simulations 
we find larger values, B — 1.9 ± 0.2 (1.3 ± 0.2) at Ki (k 2 )- 

4 Conclusions 

In this paper we have evaluated the matrix elements which contain the non- 
pertubative QCD effects in the spectator contribution to inclusive decays of the 
Ab baryon. Our principal results (in the MS scheme) are presented in table |[ The 
results indicate that spectator effects are important, accounting for a significant 
fraction of the discrepency between the theoretical prediction for r(Ab)/T(Bd) in 
eq. (0) and the experimental result in eq. ([!]) . It also appears that not all of the 
discrepency can be accounted for by spectator effects. 

The calculation described in this paper is the first evaluation, using lattice 
simulations, of the matrix elements in eq. (§) and eq. (|j) between |A 6 ) states. 
Having established that spectator effects are significant, it is now necessary to 
improve the precision, both statistical and systematic. This requires a high- 
statistics simulation at a smaller value of the lattice spacing (to decrease the 
errors due to discretisation) and with more values of the light quark mass (to 
enable a reliable extrapolation to the chiral limit). 

In this study we have used static 6-quarks. It would also be valuable, as 
a control of the systematic errors, to repeat the calculation with propagating 
6-quarks, for which these uncertainties are different. 
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